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Abstract Gronwall's function G is defined for n > 1 by Gin) = — 

•' ^ ' n log log n 

where cr(n) is the sum of the divisors of n. We call an integer > 1 
a GAl number if is composite and G{N) > G{N/p) for all prime fac- 
tors p of N. We say that is a GA2 number if G{N) > G{aN) for all 
multiples aN of A^. In arXiv 1110.5078, we used Robin's and Gronwall's 
theorems on G to prove that the Riemann Hypothesis (RH) is true if and 
only if 4 is the only number that is both GAl and GA2. Here, we study 
GAl numbers and GA2 numbers separately. We compare them with su- 
perabundant (SA) and colossally abundant (CA) numbers (first studied by 
Ramanujan). We give algorithms for computing GAl numbers; the smallest 
one with more than two prime factors is 183783600, while the smallest odd 
one is 1058462574572984015114271643676625. We find nineteen GA2 num- 
bers < 5040, and prove that a GA2 number N > 5040 exists if and only if 
RH is false, in which case N is even and > 10^^''^. 
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1 Introduction 



The sum-of-divisors function a is defined by 

d\n 

For example, a (4) = 7. 

In 1913, Gronwall [7j found tlie maximal order of a. 

Theorem 1 (Gronwall) The function 

G{n) := {n > 1) 

n log log n 

satisfies 

lim sup G{n) = e^ = 1.78107 . . . , 

n— >oo 

where 7 is t/ie Euler-Mascheroni constant. 

In 1915, Ramanujan proved an asymptotic inequality for Gronwall's func- 
tion G, assuming the Riemann Hypothesis (RH). Ramanujan's result was 
shown in the second part of his thesis. The first part was published in 1915 
while the second part was not published until much later, in 1997 



Theorem 2 (Ramanujan) // the Riemann Hypothesis is true, then 

Gin) < > 1). 

Here, n ^ 1 means for all sufficiently large n. 

In 1984, without being aware of Ramanujan's theorem, Robin [T3] proved 
that a stronger statement about the function G is equivalent to RH. 

Theorem 3 (Robin) The Riemann Hypothesis is true if and only if 

(1) G{n) < {n> 5040). 

The condition ([1]) is called Robin's inequality. Table [1] gives the twenty- 
six known numbers r for which the reverse inequality G{r) > e'^ holds (see 
[TTf Sequence A067698]), together with the value of G{r) (truncated). (The 
"a(r)" column is explained in §11 and the "Q(r)" column in §7.11 ) 

In ^14j Robin also proved, unconditionally, that 

(2) G{n) <e''+ °-^f ^ ' ^ , (n > 1) 

(log logn)^ 
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r 


SA CA GAl GA2 Factorization 


a{r)/r 


G(r) 


a(r) 


Q{r) 


3 




/ 


3 


1.333 


14.177 







4 


/ / 


/ 


22 


1.750 


5.357 





-0.763 


5 




/ 


5 


1.200 


2.521 







6 


/ / 


/ 


2-3 


2.000 


3.429 





4.134 


8 




/ 


23 


1.875 


2.561 





2.091 


9 






32 


1.444 


1.834 


4 


7.726 


10 




/ 


2-5 


1.800 


2.158 





1.168 


12 


/ / 


/ 


22.3 


2.333 


2.563 





2.090 


16 






2^ 


1.937 


1.899 


3 


1.348 


18 




/ 


2-32 


2.166 


2.041 





1.679 


20 






22.5 


2.100 


1.913 


3 


2.799 


24 


/ 


/ 


23 ■ 3 


2.500 


2.162 





1.185 


30 






2-3-5 


2.400 


1.960 


2 


1.749 


36 


/ 


/ 


22.32 


2.527 


1.980 





1.294 


48 


/ 


/ 


2^-3 


2.583 


1.908 





1.132 


60 


/ / 


/ 


22 . 3 • 5 


2.800 


1.986 





1.290 


72 




/ 


23.32 


2.708 


1.863 





1.160 


84 






22 . 3 ■ 7 


2.666 


1.791 


10 


1.430 


120 


/ / 


/ 


23 • 3 • 5 


3.000 


1.915 





1.128 


180 


/ 


/ 


22 . 32 • 5 


3.033 


1.841 





1.078 


240 


/ 


/ 


2^ • 3 • 5 


3.100 


1.822 





1.051 


360 


/ / 


/ 


23 . 32 . 5 


3.250 


1.833 





1.044 


720 


/ 




2^ • 32 • 5 


3.358 


1.782 


7 


1.028 


840 


/ 




23 • 3 • 5 • 7 


3.428 


1.797 


3 


1.065 


2520 


/ / 


/ 


23 . 32 • 5 • 7 


3.714 


1.804 





1.015 


5010 


/ / 


/ 


21-3- -S- 7 


3.838 


1.790 





1.007 



Table 1 - The set 7^ = {r < 5040 : G{r) > = 1.781 . . . }, which contains 
the subset 7^2 = {A^ < 5040 : A^ is GA2}. 
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with equality for n = 12. This refines the inequahty limsnp^^^ G{n) < 
from Gronwall's theorem. 

Recently, the authors [3] used Robin's results to derive another reformu- 
lation of RH. Before recalling its statement, we give three definitions and an 
example. 

A positive integer is a GAl number if is composite and the inequality 

G{N) > G{N/p) 

holds for all prime factors p of A^. The first few GAl numbers are 

AT = 4, 14, 22, 26, 34, 38, 46, 58, 62, 74, 82, 86, 94, 106, 118, 122, 134, 142, . . . 
(see [TTj Sequence A197638]), and (see §5.3p the smallest odd GAl number is 
N = 1058462574572984015114271643676625. 

An integer A^ > 1 is a GA2 number if 

G{N) > G{aN) 

for all multiples aN of A^. The nineteen known GA2 numbers (see Theorem [S] 
and jlTl Sequence A197369]) are 

A^ = 3, 4, 5, 6, 8, 10, 12, 18, 24, 36, 48, 60, 72, 120, 180, 240, 360, 2520, 5040. 

Every GA2 number > 5 is even. {Proof. If A^ is odd, then a{2N) = 3cr{N), 
and if A^ is also GA2, we get 

3 3G{N) _ log log 2 AT 
2 - 2G(2A0 ~ log log AT 

which implies N < 7.) 

Finally, a composite number is extraordinary if it is both GAl and GA2. 

For example, the smallest extraordinary number is 4. To see this, we first 
compute G(4) = 5.357. . . . Then, as G{2) < 0, it follows that 4 is a GAl 
number. Since Robin's unconditional bound ([2]) implies 

<^(^) < + n = 4-643 . . . < G(4) (n > 5), 

(log log 5 j"' 

we get that 4 is also GA2. Thus 4 is an extraordinary number. 

We can now recall our results from |3l Theorem 6 and Corollary 8]. 
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Theorem 4 (Caveney-Nicolas-Sondow) (i). The Riemann Hypothesis is 
true if and only if 4 is the only extraordinary number. 

(ii). If there is any counterexample to Robin's inequality, then the maximum 
fi := max{G(n) : n > 5040} exists and the least number N > 5040 with 
G{N) = fi is extraordinary. 

If there exists an extraordinary number > 4, then is even (as 5 is not 
GAl, and no GA2 number > 5 is odd) and A^ > 10*^''^ (since no GAl number 
hes in the interval [5,5040], and no GA2 number hes in [5041,10^^^^] — see 
Corollary [1]) . 

In the present paper, we study GAl numbers and GA2 numbers separa- 
tely 

Preliminary facts about GAl numbers and GA2 numbers were given in 
[3]. We recall two of them and make a definition. 

Fact 1 (proved by elementary methods in [3|, §5]). The GAl numbers with 
exactly two [not necessarily distinct) prime factors are precisely 4 and 2p, 
for primes p > 7. 

We call such GAl numbers improper, while GAl numbers with at least 
three (not necessarily distinct) prime factors will be called proper. 

The smallest proper GAl number is v := 183783600 (see §5.31 and [TTf 
Sequence A201557]). The number v was mentioned in [3j equation (3)] as 
an example of a (proper) GAl number that is not a GA2 number (because 
G{v) < G{19i^)). 

Fact 2 (see |i3i Lemma 10]). If Uq is a positive integer, then 

limsupG'(ano) = e"^, 

a—^oo 

which yields the implication 

(3) N IS GA2 =^ G{N) > e\ 

An application is an alternate proof that any GA2 number A^ > 5 is even. 
Namely, as 7 and 9 are not GA2, and as Theorem 2 in [4J says that an integer 
n > 9 is even if G{n) > e'', the result follows from 

By the method of [31 §5], one can prove two additional properties of GAl 
numbers. 

Fact 3. The only prime power GAl number N = p^ is N = 4. 

Fact 4. A product of three distinct primes pip2Pz cannot be a GAl number. 
(See §6.21 for a more general result proved by other methods.) 
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The rest of the paper is organized as follows. The next subsection esta- 
blishes notation. In §2]we recall the definitions of superabundant (SA) and co- 
lossally abundant (CA) numbers and review some of their properties. In ^we 
prove six lemmas needed later. In ^we give an analog of Theorem H] for GA2 
numbers ; in particular, if RH is false, then infinitely many GA2 numbers 
exist, and any number N > 5040 for which G{N) = max{G(n) : n > 5040} 
is both GA2 and CA. In the final four sections we study proper GAl num- 
bers : ^compares them with SA and CA numbers, ^is concerned with their 
prime factors, Ogives algorithms for computing them, and ^estimates the 
number of them up to x. 

1.1 Notation 

We let p always denote a prime. 

Let Vp{n) denote the exponent on p in the prime factorization 

n = JJp^'fW. 
p 

For n > 1, we denote the number of prime factors of n counted with 
multiplicity by 

VL{n) := Vp{n). 
p 

For n > 1, we denote the largest prime factor of n by 

P{n) := max{p : p \ n} = ma.x{p : Vp{n) > 0}. 
As usual, Chebychev's function is defined as 

p<x 

2 Review of properties of SA and CA numbers 

Superabundant and colossally abundant numbers were first introduced 
by Ramanujan, who called them generalized highly composite and generali- 
zed super highly composite numbers, respectively (cf. [1^, §59]). They were 
rediscovered later by Alaoglu and Erdos pLj. 

A superabundant (SA) number is a positive integer N such that 

a(N) a(n) 

——^ > (0<n< N). 
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The first few SA numbers are (see [T7t Sequence A004394]) 

= 1, 2, 4, 6, 12, 24, 36, 48, 60, 120, 180, 240, 360, 720, 840, 1260, 1680, .... 

A colossally abundant {CA) number is a positive integer A^ for wliich 
tliere exists an exponent £ > sucli that 

(4) . 

Such an exponent e is called a parameter of A^. The sequence of CA numbers 
(compare [T71 Sequence A004490]) begins 

A^ = 1, 2, 6, 12, 60, 120, 360, 2520, 5040, 55440, 720720, 1441440, 4324320, .... 

From (jlj), it is easy to show that every CA number is also SA. 

Now let A^ denote an SA or CA number. Then (see jU Theorems 1 and 3] 
or [la §59]) 

(5) AT = 2^=2 . 3^3 . . . . pfcp ^ k2>h>h>---> kp 
with kp = 1 unless A^ = 4 or 36, and [1, Theorem 7] 

(6) p = P(A^) ~ logA^ (A^^oo). 

We recall some properties of CA numbers (see [H [21 El El El [ISl [HI [U])- 
Note first that for any fixed positive integer k, the quantity 



log ( 1 + ^ 



F{t,k) :-- 



t + f^ + --- + t^ 



logt 



is decreasing on the interval 1 < t < oo, and the function t i— )■ F(t, k) maps 
the interval onto the positive real numbers. Hence, given £ > 0, we may 
define Xk = Xk{e) > 1 by 

(7) F{xk,k) = e. 

(See [T^ p. 189] and [T3| §61 and §69].) In particular, when k = 1 we set 
X = xi = xi{e), so that 



(8) F(x,l) = F(xi,l) 



log(l + - 

\ X 

logx 



It is convenient to set Xq = +00. From the decreasingness of F{t, k) with 
respect to both t and k, it follows that the sequence {xk)k>o is decreasing. 

If is a CA number of parameter e and p divides with Vp{N) = k, 
then applying (jlj) with n = Np yields 

e > F{p, k + 1) i.e. p > Xk+i 

while, if > 0, applying with n = N/p yields 

(9) e < F{p, k) i.e. p < xt- 

Let K be the largest integer such that xk > 2. Then from ([9]), for all p's we 
have 2 < p < Xk and 

k = Vp{N) < K. 

Now define the set 

S := {F{p, k) : p is prime and k > 1}. 
Its largest element is 

max £ = F(2, 1) = = 0.5849 . . . , 

log 2 



and its infimum is 



inf S = lim F{p, k) = 

fc— >oo 



for any fixed prime p. 

If e ^ S, then no Xk is a prime number and there exists a unique CA 
number A^ = N{e) of parameter e; moreover, A^ is given by either of the 
equivalent formulas 

AT = JJ with Xk^+i <p<Xk^ 

p<x 

or 

K 

(10) N = l[ Hp. 

k=l p<Xk 

In particular, if e > max£^, then x = xi < 2, K = and N{e) = 1. 

U e & S, then some Xk is prime, and it is highly probable that only one Xk 
is prime. But (see ^ Proposition 4]), from the theorem of six exponentials 
it is only possible to show that at most two x^s are prime. (Compare ^ 
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p. 538].) Therefore there are either two or four CA numbers of parameter e, 
defined by 



K 



(11) 



k=l p<X], 



or 



p<Xk 



Here, if Xk is a prime p for some k, then p may or may not be a factor in the 
inner product. (This can occur for at most two values of k.) In other words, 
if Xk-i < p < Xk, then the exponent Vp{N) of p in N is k, while if p = Xk, the 
exponent may he k ot k — 1. In particular, if is the largest CA number of 
parameter e, then 



Note that, since if e ^ S, then Xk is not prime, formula (fTOj) gives the same 
value as f llip . Therefore, for any e, formula f llip gives all the possible values 
of a CA number A^ of parameter e. (Thus A^ is a product of "primorials" \T7\ 
Sequence A002110].) 

3 Six lemmas 

The case k = 2 of the following lemma was proved in |14[ p. 190]. 
Lemma 1 For k > 2, we have the upper bound 



Proof. Since the function t F{t, k) is strictly decreasing on 1 < t < oo, to 
prove Xk < z := [kxY^^, it suffices to show F(z,k) < F(xk,k). As ([7]) and 
([H]) imply F{xk, k) = e = F{x, 1), this reduces to showing F{z, k) < F{x, 1). 
Since z > 1 and k >2, we have 



(12) 





Xk < {kxY^^. 




1 k 



< 



{z + z'^ + ■ ■ ■ + z'') log z {z + z'^ + ■ ■ ■ + z^) log kx 



k k 1 



< 



{k - 1 + z'')\ogx (| + /cx)logx (x + |)logx 



< 
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using the lower bound log (l + j) > (t + |) ^ , valid for t > 0. This proves 
the desired inequality. □ 

In the proof of Theorem [S] (iii), we will need the following result (see [TTj 
Lemma 4]). 

Lemma 2 Given a CA number Nq of parameter Eq, let N > Nq be a number 
satisfying 

, , a(n) (t(N) 

13 n> No ^ ^ < 

for some fixed e > 0. Then N is CA of parameter e. 



Proof. Since Nq is CA of parameter e^, we have 

l+eo 



On the other hand, (fT^ yields 



a{N) ^ f N 



a(iVo) -\No 



l+e 



Hence e < Eq. 

In view of f fT3|) . to prove that is CA of parameter e, we only need to 
show that 

If n < Aq, then since A'o is CA and (fT3|) holds, we have 

a{n) _ ^ a{No)n''''' ^ (t(A'o) A'q""^ _ a(A'o) ^ f^(A^) 

This completes the proof of Lemma |2j □ 

The next lemma provides an estimate for a CA number of parameter e. 

Lemma 3 Let N be a CA number of parameter e < F(2, 1) = log(3/2)/ log2 
and define x = x{e) by ([8]). 
(i). Then 

log A^ < (x) + cy/x 

for some constant c > 0. 
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(ii). Moreover, if N is the largest CA numhei^ of parameter e , then 

(x) < log N < {x) + cy/x. 

Proof, (i). It follows from formula (fTT]) for that if Xk is defined by ([7]), then 

(14) \ogN <{xi) + {x2) + ■■■ + {xk), 

where K is the largest integer such that xk > 2. (Note that V2{N) = K or 
K - I, and that e < F(2, 1) implies x > 2 and K > 1.) 
As t ^-)■ F{t, k) is decreasing and ([7]) holds, we have 

F{2, K) > F{xK, K) = e = F{xk+u K + 1) > F{2, K + 1). 

On the other hand, 

F{2,K) = log (^1 + 2X+i_2) ^ (2^+1 -2) log 2 " 2^n^ ^ 2^ 

and, from ([8]), 

log(l + i) 1^1 1 

e = — ^ — > > > — . 

logx (x + l)logx (x + l)(x— 1) 

Thus 

A>F(2,ir)>.>l, 

implying 

2 

(15) K<1 + - logx. 

log 2 

Since k ^ x^ is decreasing, from fll4p we have (compare [T3, equation 
(368)]) 

log AT < (Xi) + (X2) +ir(x3). 

Using X2 < a/2x and X3 < \/3x (from Lemma [1]), together with (fTSl) and the 
Prime Number Theorem in the form (t) ~ t, we deduce (i). 

(ii). From (fTTl) . the largest CA number of parameter e is 

k=l p<Xk 

1. Note that Ramanujan's definition of CA number of parameter e in |13j is not exactly 
the same as that of Robin in [14, pp. 189-190]. Ramanujan's definition corresponds to the 
largest CA number of parameter e for Robin. 
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which imphes (x) < logA^, and (ii) follows from (i). □ 

In the next lemma, we recall the oscillations of Chebychev's function 
studied by Littlewood. 

Lemma 4 There exists a constant c > such that for infinitely many 
primes p we have 

(16) 9{p) < p — c^/p\og\og\ogp, 
and for infinitely many other primes p we have 

(17) 9{p) > p + c^/p\og\og\ogp. 

Proof. From Littlewood's theorem (see [lOj), we know that there exists a 
constant c' > such that for a sequence of values of x going to infinity we 
have 

(18) 9{x) < X — c'A/xlogloglogx, 

and for a sequence of values of x' going to infinity we have 

(19) x' + CAi/a/logloglogx' < 6'(a;'). 

Let us suppose first that x is large enough and satisfies (llSp . If x = p is 
prime, then (1181) implies (1161) . Now assume x is not prime, and let p be the 
prime following x. As the function t ^ t — cy^ log log log t is increasing, we 
get 

0{j)) = 9{x) + \ogp < X — c' log log log X + logp 
< p — c'-y/plogloglogp + logp, 

which implies (I16p with c < c' for x large enough. 

The proof of (11 7p is easier. Let x' satisfy (I19p and choose the largest prime 
p < x' . For c < c', we have 



9{p) = 9 (x') > x' + c V x' log log log x' > p + Cy/p log log log p, 
which proves (fT7|l . □ 

Lemma 5 Chebychev's function (x) satisfies 

(x) < (1 + a)x, 

where 

z/x<8-10i\ 

CX = Qi(x) '= \ 1 

^ < 0.000028 otherwise. 

36260 
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Proof. Schoenfeld (cf. [TB| p. 360]) proved {x) < 1.000081a; for all x, and he 
mentioned that Brent had checked that (x) < x for x < 10^^. The stronger 
results stated here are due to Dusart — see [5l p. 2 and Table 6.6]. □ 

Lemma 6 Let e be a positive real number. For t > e, let us set 

(20) g{t) = g,{t) := elogt - log log log t. 

Then there exists a unique real number t^ = tole) > e such that 

(21) ^ 



log to log log to 

Moreover, g{t) is decreasing for e < t < to and increasing for t > t^. 
Proof. The derivative of g is 

g'it) = Ue- ^ 



t \ logt log log t 

For t > e, both log t and log log t are positive and increasing, and the function 
t I— i- l/(logt log log t) is a decreasing bijection from (e, +oo) onto (0, +oo). 
Therefore, one can define to > e by fl2ip . 

Then we have g'{t) < for e < t < to? and g'{t) > for t > to? which 
completes the proof of Lemma [6l □ 



4 GA2 numbers 

We first study GA2 numbers. Compare the following result on them with 
Theorem H] on extraordinary numbers. 

Theorem 5 (i). The set of GA2 numbers < 5040 is 

7^2 := {3, 4, 5, 6, 8, 10, 12, 18, 24, 36, 48, 60, 72, 120, 180, 240, 360, 2520, 5040}. 

(ii) . If the Riemann Hypothesis is true, then no GA2 number exceeds 5040. 

(iii) . If the Riemann Hypothesis is false, then infinitely many GA2 numbers 
exist; moreover, the inequality 

/i := max{G'(n) : n > 5040} > 

holds, and any integer A > 5040 for which G{A) = fj, is both GA2 and GA. 
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Proof, (i). Setting 

7^' := {N < 5040 : is GA2}, 

we have to prove that TZ' = 7^2- 

To show TZ' C 7^.2, choose N G TZ'. From ([3]), we have G{N) > e'^ , so that 

N en:={r< 5040 : G(r) > e^} 

= {3, 4, 5, 6, 8, 9, 10, 12, 16, 18, 20, 24, 30, 36, 48, 60, 72, 84, 120, 180, 
240, 360, 720, 840, 2520, 5040}, 

by calculating the "r" column of Table [T] To show that belongs to the 
subset 7^2 C TZ, define for r G 7^ the integer 



a(r) := 



min Ar, if Ar := {a : G{ar) > G{r), ar G TZ] ^ 
0, if Ar = 0. 



A computation (see the "a(r)" column of Tabled]) shows that 

(22) {r G 7^ : A ^ 0} = {9, 16, 20, 30, 84, 720, 840}. 
Since A^ is GA2, it must lie in the complement 

(23) 7^ \ {9, 16, 20, 30, 84, 720, 840} = 7^2. 

This shows TZ' C 7^2- 

To prove 7^2 C TZ', choose r G 7^2- To get r E TZ', we need to show that 
G{r) > G{ar), for any multiple ar of r. We consider two cases. 

Case 1 : ar < 5040. If ar G TZ, then since r G 7^.2, relations (123|) and f l22|) 
imply G'(ar) < G{r). On the other hand, if ar ^ 7^, then G{ar) < e"* < G{r). 
Thus G{r) > G{ar) whenever ar < 5040. 

Before considering Case 2, we recall that in [T^t p. 204 (c)] Robin proved 
that if C is the largest CA number with P{C) < 20000, then there is no 
counterexample < C to his inequality ([T]). From the property of CA 
numbers, we have logC > (20000), where (x) is Chebychev's function. 

We also recall that in |16[ p. 359, Corollary 2], Schoenfeld proved that 

(x)>x-—^ (x> 19421). 

A calculation then gives the inequalities 

(20000) > 20000 > 19747 

^ ' 8 log 20000 
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which, together with Robin's resuh on C, yield the imphcation 
(24) 5040 <n< e^^^^^ =^ G{n) < e^. 



Case 2: ar > 5040. If log ar < 19747, then ^ gives G{ar) < < G{r). 
On the other hand, if log ar > 19747, then from ([2]) we get 

G{ar) <e^+ °f„^^f-., = 1.787 . . . < 1.790 . . . = min G(r') < G{r). 
(logl9747)'' r'en2 

Thus G{r) > G{ar) whenever ar > 5040. 

This shows that, in both Cases 1 and 2, all elements r of 7^2 are GA2 
numbers, so that TZ2 C TZ'. Finally, since we already have TZ2 D TZ', we get 
7^2 = TZ'. This proves (i). 

(ii) . If RH holds, then by Robin's theorem there is no number n > 5040 with 
G{n) > e"', while from a GA2 number must satisfy G{N) > . 

(iii) . Let us assume that RH fails. Set 

e := sup 5R(p) 
C{p)=o 

so that 

1/2 < e < 1. 

Let denote a CA number of parameter e, and define x = x{6) by (jSj). 
If p := P{N) and if p"*" is the prime following p, then from (lll|) we have 

p < Xi = X < , 

which implies x ~ p as A^ — > oo. Further, from we get p ~ log A^, which 
implies 

X ~ log A^ {N ^ oo). 
In [T5l p. 241], it is proved that as A^ — )■ oo 

G{N) = e'' {1 + Q+ {x-'')) (l-e<6<l/2) 

which implies that 

G{N) = (1 + Q+ ((log AT)-'')) (1 - 6 < 6 < 1/2). 

(Here the notation "f{N) = n+{g{N)) as N oo" means that f{N) > g{N) 
infinitely often, and should not be confused with the notation Q{n) in §1.11 ) 
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Therefore, there exist infinitely many CA numbers satisfying G{N) > e'^, 
and, for all t, we have ma.Xn>tG{n) > . 

Now we construct two sequences Ai, ^2, . . . and A[, A'2, . . . , as follows. 
Let Ai (resp., A[) be the smallest (resp., largest)[l integer > 5040 such that 
G{Ai) = gIa[) = fi. 

Given i > 2, assume that Ai, A2, . . . , Ai_i and A[, A2, . . . , A^_i have 
been defined. Set /ij := max„>yi/ G{n) and let A^ (resp., A'^ be the smallest 
(resp., largest) integer > A^_i with G{Ai) = ^(^4-) = /i,. Since we have 
fii > e"' = lim sup (?(«), infinitely many Aj's can be found. The numbers Ai 
are such that 

n> Ai ^ G{n) < G{Ai) 

and, therefore, are GA2. 

In the same way, A is proved to be GA2, using A > 5040 and G{A) = /i. 
To show that A is CA, we apply Lemma [2] with Nq = 55040, Eq = 0.03, 
N = A, and e = l/(logy41oglogy4) ; since A is GA2 and A > 5040, from 
fl2i|) and (jS]) we obtain that N = A > e^^^^^ > A'o- For > A'o, from the 
definition of A we have G{n) < G{A). Since e < Aq < ^ holds, it follows 
from Lemma [H] that, on the interval [A"o, +C)o), the function g{t) (defined by 
(120|) ) attains its minimum at t = A. Thus, for n > Nq, we have 

4^ = G{n)^-^^ = G(n)e-W < G(A)e-('^) < G(A)e-(^) = ^ 

and so f|T3|) holds. Applying Lemma |2] completes the proof of (iii). □ 
Here is a corollary of the proof of Theorem [5l 

Corollary 1 There is no GA2 or extraordinary number between 5041 and 

]^g8576 

Proof. Since lO^^^^ < 6^9747^ ^^As follows from □ 

5 Comparison between CA and GAl numbers 

In this section, we study GAl numbers. We begin by comparing them 
with CA numbers. 

2. It is highly probable that Ai — A[. A difhcult question is whether G is injective. 
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5.1 CA and GAl 

By revisiting the proof of [T5| Theorem 3, p. 242], we shall prove the 
following results. 

Lemma 7 Let N be a CA number of parameter e > and assume that 
P ■■= PiN) > 5. // 

(25) e > ^ 



\og{N/p) log log{N/p) ' 
then N is also a GAl number. 

Proof. Let g be a prime factor of N. It follows from (|5]) that 6p divides N 
and that N/q > N/p > 6 > e, which implies loglog(A^/g) > log log e = 0. 
Since is a CA number, from (jl]) one has 

^ ^ {N/qy+^ - 

so that 

(27) ^< ' 



Since log log and loglog(A^/g) are positive, it follows that 

^ qnoglogiN/q) = ArMoglog(Ar/g) = ^-^(^iN/q) - g{N)), 

where g{t) is defined by (l20l) . By Lemma [6], using (12T]) to define to > e, we 
have that g{t) is increasing for t > to. Now from (125|) we deduce that 

e<to< — < — <N 
p q 

and from 1^ we get G(A^/g) < G{N). This shows that A^ is GAL □ 
Theorem 6 Infinitely many CA numbers are GAl. 

Proof. Choose a sufficiently large prime p satisfying f lT7|) . and set e := F{p, 1) 
(so that X = p, hj (^). Let A^ be the largest CA number of parameter e (so 
that p divides A^, by (fT2|) ). From Lemma [3] part (ii) and (fT7|) . we get 

log N > {x) = (p) > p + cy^log log logp. 
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so that 

log(A^/p) > p + c^/p log log log p — logp > p + 1. 
Using the lower bound log(l + t) >t/{l + t), we get 



1 



log (l + i 
Fip, 1) = — ^ ^ > 



logp {p + 1) logp 



> : r > 



{p + 1) log(p + 1) log(iV/p) loglog(iV/p) 

and Lemma [7] implies is GAl. Since, by Lemma HI there are infinitely 
many primes p satisfying (fT7|) . the theorem is proved. □ 



5.2 CA and not GAl 

To study CA numbers that are not GAl, we need a lemma. 

Lemma 8 Given a prime p > 3, let N be the largest CA number of parame- 
ter e := F{p,l). If 

(29) e < ^ 



log A^ log log A^ 
then N is not GAl. 

Proof. As e = F{p,l), we have p"^ = {p + l)/p = cr{p)/p- Hence, by (fT2l) . 
inequality (126|) becomes an equality when q = p, and so do inequalities ( 127|) 
and ( !28|) . Therefore, with g and to defined by (120|) and (12T|) as in the proof 
of Lemma [TJ we get that 

^^^ = exp{g{N/p)-g{N)) 

and, from Lemma El that g{t) is decreasing for t < t^. Then (12^ implies 
N/p <N <to, so that G{N) < G{N/p). Thus A^ is not GAl. □ 

The CA numbers A^ such that P{N) e {2,3,5,7,11,13,29,59,149} are 
not GAl. There are two CA numbers such that P{N) = 23 ; the larger one 
is not GAl, while the smaller one is GAl. All other CA numbers satisfying 
P{N) < 300 are GAl. (These statements follow by computing all CA num- 
bers A^ with P{N) < 300, and calculating those that are GAl — see §7.) 

Theorem 7 Infinitely many GA numbers are not GAl. 



18 



Proof. Choose a sufficiently large prime p satisfying f lTO]) . and set e := F{p, 1) 
(so that, from ([H]), x = p). Let be the largest CA number of parameter e 
(so that, from (fT2|) . p = P{N)). From Lemma [3] part (i) and (fT6l) . we get 

log N < (p) + < p — cy^ log log logp + Ci/p < Pi 

and so 

log(l + ^) 1 1 

^ - < ^ < 



logp plogp log log log 

Then Lemma [8] implies A^ is not GAL Since there are infinitely many primes p 
satisfying (fT6|) . the theorem is proved. □ 



5.3 Odd GAl numbers 

We show that there are infinitely many odd GAl numbers, and we com- 
pute the smallest one. 

Let us denote by Vq = {2, 3, 5, 7, 11, 13, 17, ... } the set of all primes, and 
by P a subset of Pq- To P, we attach the set 

A/p := {n > 1 : p \ n =^ p G V} 

and the function 

:= Yl 

PdiP, p<X 

A number A^ G N'-p is said to be colossally abundant relative to V (for 
short, CA-p) if there exists e > such that 

a{N) ain) 

If M = ripgPo ^"'^ ordinary CA number of parameter e, then the factor 
A^ = ripeP^'"'' CAp, for the same parameter e, and all CA-p numbers can 
be obtained in this way. 

Theorem 8 There exist infinitely many odd GAl numbers. 

Proof. First, we observe that Lemma [7] remains valid if we replace CA with 
CAp, for any set V with at least 2 elements. 

We set V = Vo\ {2}. The proof of Theorem [6] remains essentially valid. 
We just have to change the lower bound for log A^ to 

logN >p{p) = {p)-\og2 
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and the inequality log(A^/p) > p+1 still holds, so that we may conclude that 



The smallest CApg\{2} number that is GAl is 

CO := 1058462574572984015114271643676625 

= 3^-5^-7^-ll^-13-17-19-23-29-31-37-41-43-47-53-59-61-67-71-73. 

From our computation (see §7.5p . u is also the smallest odd GAl number. 

Corollary 2 There exist infinitely many GAl numbers that are not SA. 

Proof. This folllows immediately from ([5]) and Theorem [HI □ 

Of course, the proof of Theorem |8] works for any set of primes V such 
that Vo\V is finite. 

6 Prime factors of GAl numbers 

Here we study prime factors of proper GAl numbers. 

6.1 An upper bound 

We need the following upper bound. 

Theorem 9 Given a GAl number N with Q{N) > 3, letp be a prime factor 
of N. Then for any positive integer r < Vp{N) we have 



N is GAl. 



□ 



p < {r log A^) 



< logA^. 



Proof. We have G{N/p) < G{N), which implies 



(30) 



a{N/p)N 
{N/p)a{N) 



< 



loglog(A^/p) log (log A^ — log p) 
log log A^ log log A^ 



Note that log log A^ > loglog(A^/p) > log log 4 > 0. We also have 



log (log A^ — logp) 
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so that 

log(log iV - logp) _ ^ ~ ~ toff^) 

log log log log N 

Further, setting v = Vp{N), the left side of fl30p can be written as 



a{N/p)N Vp 



v-l 



{N/p)a{N) ^ 1 + p^ hp" 

1 1 

32 =1 > 1 

From dsn]), (EU), and ([22]), one deduces 

(33) f<l+p + ...+f< , ^ log AT log log iV 

-log(l-S) 

which yields 

(34) p^ log p < log log log A^. 

Let us assume, ah absurdum, that p > (r log NY^^ . Then we would have 

p^ logp > (r log A^) - log(r log A^) = log A^ log(r log A^) > log A^ log log A^ 

r 

contradicting (IMj) . Therefore, p < {r log NY^^ holds. Finally, by calculus, 
(r log NY^^ is decreasing for r > 1 (because n{N) > 3 implies A^ > 8 and 
log N > 2) and the theorem follows. □ 

6.2 Study of n{N) where N is GAl 

We show that there are only finitely many proper GAl numbers A^ that 
have a fixed value of Q{N). 

Theorem 10 If k > 3, then 

Uk := #{A^ : A^ is GAl and n{N) = k} < oo. 

Proof. For a GAl number A^ with Q{N) = k > 2, let us write A^ = pip2 ■ ■ - pk 
with pi < p2 < ■ ■ ■ < Pfc- We have A^ < P^j so that pk > N^/'^ holds. But 
Theorem [9] yields pk < log N, whence 



log log N 
21 



and is bounded. Thus is finite. □ 

Since J°f„g°o6o = 28.03 . . . , a table of GAl numbers up to 10^° (see §7]) 
allows us to calculate 11^ for k < 28. 

We have 11^ = if 3 < A; < 12, and the following table gives 11^ when 
13 < A; < 28 (see [17,, Sequence A201558]). 



k 



13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 
2421124123777147 



6.3 The exponent of the largest prime factor 

First, we observe that the function t ^ 2*/t is an increasing bijection of 
the interval [2, +oo) to itself. Let us introduce the inverse function h defined 
for X > 2 by 

2* 

(35) h{x) = t <^ J' 
We shall need the following lemma. 

Lemma 9 Let x satisfy x >2. Then we have 2 < h{x) < 3.08 logx. 

Proof. The lower bound results from the definition of h. Let us set t = h{x), 
so that X = 2^/t. By noting that (logt)/t < 1/e holds, we get 

h(x) t 1 1 
^ ^ = = < = 3.0743 

logx tlog2 — logt log2 — (logt)/t ~ log2 — 1/e 
which proves Lemma O □ 

Theorem 11 Let N be a GAl number with Q{N) > 3. Set R = h{\ogN), 
so that = logiV. Then N divides the number M = M{N) defined by 

[R] LRJ 

(36) M:=n n P=UUp 

^•=1 ((T-+l)logAf)l/('-+l)<p<(rlogAf)i/'- r=l p<(r log Af)!/'- 

Proof. Since the function r y-^ (r\og N)^^"^ is decreasing, this follows from 
Theorem O □ 

For example, ii N = u = 183783600, we compute R = h{log v) = 7.072 . . . 
and find that M = 72u. 



22 



Theorem [TT] allows the computation of proper GAl numbers — see §7.2 
and §7.5. 

For the exponent Vp{N) of a prime p in the standard factorization of N, 
Theorem [TTl provides the upper bound Vp{N) < Vp{M), which only depends 
on the size of A^. 

We now study the exponent of the largest prime factor of a GAl number. 

Theorem 12 Let N be a GAl number with ^1{N) > 3, and let p = P{N) be 
its largest prime factor. Then Vp{N) = 1. 

Proof. Suppose on the contrary that v := Vp{N) > 2. Then Theorem [9] implies 
that A^ divides the number 



M„ = M,{N) := I n 

yp<{v logAf)!/" 



n n p 

r=v+l p<{r log Af)!/'- 



with R defined by 2 /R = \ogN. Thus, from the function r i— (rlogA^) 
being decreasing, 



l/r 



logN < logM^ = v{v\ogN)^/'')+ J2 {{rlogNy/') 

r=v+l 

< 2{{2\ogNy/^) + R{{3\ogNy/^). 
From Lemmas [H] and it follows that 

logN < 1.000028 (2^/2 log A^ + 3.08 log log N{3 log N)^/^ 
< 2.83v/logiV + 4.45 log log A^(logA^)^/^ 
Therefore, we have 

2.83 4.45 
y/hiW ^ (log Ar)2/3 log log A^ - ^ 

which implies logA^ < 15.03, A^ < A^o := 3336369 and R < 6.65, so that 
A^ must divide M„(A"o) for some v in the range 2 < i; < 6. But the table 



V = 


2 


3 


4 


5 


6 


{vhgNoy/'" = 


5.48 


3.56 


2.78 


2.37 


2.12 


M^iNo) = 


43200 = 2^3352 


1728 = 2633 


64 = 2^ 


64 


64 



shows that if w > 2, then the number M„(A'o) divides M2(A'o) = 43200, 
contradicting the easily-checked fact that none of the 84 divisors of 43200 is 
a proper GAl number. (In fact, we will see in §7.51 that there is no proper 
GAl number < 183783600.) This proves the theorem. □ 
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6.4 The largest prime factor of a GAl number 

For a GAl number, we now study the largest prime factor itself. 

Theorem 13 For GAl numbers N with Q{N) > 3, the largest prime factor 
satisfies 

P{N)r^ log N (N^oo). 

Proof. Let be a GAl number satisfying Q{N) > 3 and let p := P{N) be 
its largest prime factor. From Theorem [9l we know that 

(37) p<logA^. 

It remains to get a lower bound for p. The proof resembles that of Theorem [T2l 
Since A^ divides M given by (136|) and p = P{N), by Lemma [5] we have 

R 

logN < (p) + ^((rlogAr)V'-) 

r=2 

< (p) + ((21ogAr)i/2) + o(i?(iogiV)i/3) 

(38) < {p) + 0{^/h^). 

From the Prime Number Theorem and from (137|) . we get 

(p) = p + 0{pexp{— c\/ log p)) = p + (9(log A^exp(— c^/log log A^)). 
Therefore, fl55]) becomes 

logA^ < p + 0(logA^exp(-CA/log log A^)), 
which, together with (1371) . completes the proof of the theorem. □ 

7 Computation of GAl numbers 

In this section we give several versions of an algorithm to compute GAl 
numbers. 

7.1 The Gronwall quotient 

We begin with a lemma and a definition. 
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Lemma 10 Let n be a positive integer with Q{n) > 3. Let q and p be prime 
factors of n satisfying q < p and Vq{n) < Vp{n). Then we have 



G{n/q) < G{n/p). 



Proof. We have 



cr{n/q)/{n/q) cr{n/q) q -\ \- qM"') i 



< 1 



a{n)/n a{n) 1 + gH ^ qM-^) 1 + gH \- qM^^) 

1 (^{^/p)/{tt-/p) 
1 + p + ■ ■ ■ + p^p^'^'> a{n)/n 



which imphes 21!]:/^ < , 

n/q — n/p 

The lemma follows from loglog(n/g) > loglog(n/p), since VL{n) > 3 im- 
plies log log (n/p) > log log 4 > 0. □ 



We define the Gronwall quotient Q{n) of a composite integer n to be the 
number 

^, , Ginip) p"p(-n)+^ — p loglogra 

Q{n) := max = max 



p\n G{n) p\n p«'p(")+i — 1 loglog(n/p) 

p prime p prime 

GAl numbers are characterized by Q{N) < 1. For example, the "Q(r)" 
column in Table [T] shows that the only GAl number r G 7?. is r = 4. 

Let us introduce a subset S{n) of the set of the prime divisors of n. The 
elements of S{n) are defined by induction. The largest prime factor of n is the 
first element qi of S{n). Now let us assume that i>2 and that the elements 
gi, q2, . . . , qi-i G S{n) are known. 

If, for all primes p that divide n and are smaller than we have 

Vp{n) < Vq^_^{n), then there are no further elements of S{n), and we get 
S{n) = {quq2,. . .,qi-i}. 

If not, then the element qi G S{n) is defined as the largest prime factor 
of n that satisfies qi < gj_i and Vq.{n) > 

From Lemma [TOl if Q{n) > 3 we get 

^, , Gin/p) pfpW+i _ j9 loglogn 

Q[n) = max = max 



pe5(n) G{n) pG5{n) p^p 1 loglog(n/p) 
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7.2 A first algorithm 

To compute all proper GAl numbers N < x for a given x, we first calcu- 
late M = M{x), defined by 

M:=n n ^ 

with R such that [Rlogxy/^ = 2. Any GAl number N <x with Q{N) > 3 
is a divisor of M (see Theorem [TT]) . 

Thus a first version of the algorithm computes all composite divisors 
of M, and for each of them calculates G{N/p)/G{N) for all p G S{N). If 
for some p e S{N) we have G{N/p)/G{N) > 1, we stop : A^ is not GAl. If 
not, we compute the Gronwall quotient Q{N) (which involves all primes p 
dividing A^) : A^ is GAl if and only if Q{N) < l.B 

7.3 A second algorithm 

A more elaborate version of the algorithm tests only a small number of 
the divisors of M. First, we define 

M,:=n n p 

r=2 p<(^r logx)^/^ 

SO that M2 := M/Mi is squarefree. Let us write M2 = PiP2'''Ps where 
Pi,P2, ■ ■ ■ ,Ps are consecutive primes in ascending order. 

As a first step, we compute the set Vq of all the composite divisors of Mi 
and test each of them for GAl by the method described above. 

A divisor of M whose largest prime factor is pi is equal to dpi, where d is 
a divisor of M whose largest prime factor is < Pi. Therefore, we construct by 
induction on 2 = 1, 2,. ..,s the set V- containing those divisors of M whose 
largest prime factor is Pi, and the set "Dj containing the divisors of M whose 
largest prime factor is < pi. Then V[ is equal to piDi^i and Vi ='D[U I^j-i. 
From Theorem [9l for z = 1, 2, . . . , s, we only have to test the elements of T>[ 
that are greater than exp(pj). 

3. To avoid roundoff errors, we carry out our computation in floating point arithmetic 
with 20 decimal digits and choose a smaU e (typically, e = 10^^). In the first step, we keep 
the A^'s satisfying Q{N) < 1 + e. For these A^'s, we start the computation again with 40 
digits. 
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7.4 A third algorithm 

Let us say that a divisor d G T>i (with < i < s) is bad if, for every 
j satisfying i < j < s, all multiples of d belonging to Vj are smaller than 

The largest multiple of d belonging to Vj is dpi^iPi^2 ' ' 'Pjj so that d is 
bad if and only if 

logo? < 6i := (pi) + min (pj - (pj)). 

i<j<s 

Therefore, we write Qi C Vi for the set obtained from by removing the 
bad divisors, i.e., those divisors d satisfying c? < A, := exp(5j). 

Furthermore, we construct Q'i^i and QiJ^i by removing from Pi+iQi and 
Pi+iQi U Qi, respectively, those divisors d that satisfy d < Aj+i = exp(5j+i). 

For z = 1, 2, . . . , s, it remains to test the elements of Qi whose largest 
prime factor is equal to pi, that is, the elements of Q[. 

7.5 Results 

The smallest proper GAl number is 

V = 183783600 = 2^ ■ 3^ ■ 5^ ■ 7 ■ 11 ■ 13 ■ 17. 

We compute that M = M(z/) = 8 ■ 19 ■ u and we find that there is no proper 
GAl number N < u. 

Using the third algorithm, we have computed all GAl numbers < 10^° 
with n{N) > 3. 

These results as well as the Maple code can be found on the web site 
http : / / math . un iv-lyonl .f r/~nicolas/GAnumbers .html, 

We hope to present soon a fourth algorithm, more sophisticated, and able 
to compute GAl numbers up to 10^^°. 

8 The number of GAl numbers up to x 

Let Qi{x) be the number of proper GAl numbers N < x. From (I36p we 
know that Qi{x) does not exceed the number r(M) of divisors of 

[Rj 

M = M{x) = W W p 

1'=^ p<{r log x)^^^ 
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with (RlogxY^^ = 2. It is easy to see that logM ~ logx as x — )■ oo, and 
from the estimation of the large values of the function r (cf. [S] or it 
follows that 

Qi(x) < exp (c "f^ ) 
\ log log X J 

for some positive c. By estimating the number of good divisors of M (that 
is, divisors that are not bad — see §7.4p . it might be possible to improve the 
above estimate. 

It seems more difficult to get a lower bound for Q\[x). We hope to return 
to these questions in another article. 

Acknowledgement The authors thank Kieren MacMillan for his help with 
TeXing the manuscript. 
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